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Monte Carlo Simulations for Top Pair and Single Top Production 
at the Tevatron 
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(J— |) Deutsches Elektronen- Synchrotron, Platanenallee 6, 15738 Zeuthen, Germany 



Summary. — Monte Carlo (MC) simulations are indispensable tools for top quark 
physics, both at the current Tevatron collider and the upcoming Large Hadron Col- 
lider. In this paper we review how the Tevatron experiments CDF and D0 utilize 
MC simulations for top quark analyses. We describe the standard MC generators 
used to simulate top quark pair and single top quark production, followed by a dis- 
cussion of methods to extract systematic uncertainties of top physics results related 
to the MC generator choice. The paper also shows the special MC requirements for 
some example top properties measurements at the Tevatron. 

PACS 14.65. Ha - Top quarks. 

PACS 2.70.Uu - Applications of Monte Carlo methods. 



1. Introduction 

With more than 3.5 fb _1 of data recorded per Tevatron experiment to date (Sum- 
mer 2008), top quark physics has entered an era of precision measurements. The latest 
combination of top quark mass measurements from the CDF and D0 collaborations 
shows a relative uncertainty of only 0.8% [1], and multivariate analysis techniques are 
being used in many analyses, notably the recent evidence for the electroweak production 
of single top quarks [2,3]. For all these endeavors, Monte Carlo (MC) simulations are 
key tools, for estimating the signal acceptance and the amount of background, and for 
evaluating systematic uncertainties. 



2. Standard Monte Carlo Generators at the Tevatron 

2T. Top Pair Production. - The Tevatron collaborations follow different approaches 
for their default MC generators for top quark pair production, as summarized in table [U 
Both collaborations use leading order (LO) MC generators with parton showering (PS), 
Pythia v6.2 [4] in the case of CDF, and ALPGEN v2.1 [5] (using Pythia v6.3 for PS) 
for D0. As both MC generators are LO generators, the generated tt production cross 
sections must be scaled to the theoretical cross section. CDF has opted for a well estab- 
lished MC generator that has been carefully validated against and tuned to the Tevatron 
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Table I. - Comparison of the default MC generators used for top quark pair production by the 
CDF and D0 collaborations. 

CDF D0 

Generator Pythia v6.2 ALPGEN v2.1 with Pythia v6.3 

Process qq, gg — * ti qq, gg — > tt+ 0-2 partons 

Parton Distribution Functions CTEQ5L [6] CTEQ6L [7] 

Tunes Tune A, W/Z p T 

Multiple Collisions Minimum Bias (Pythia) Zero Bias (Data) 



data. D0 uses a fairly recent MC generator with additional functionality compared to 
Pythia, including exact matrix elements for 1 ^ n processes with a matching procedure 
between partons generated via matrix elements and parton showers, a more recent set of 
parton distribution functions (PDFs), and it spin correlations. The collaborations have 
also chosen different ways of simulating multiple collisions in the same bunch crossing. 
CDF overlays minimum bias events generated with Pythia, with the same calibration 
constants applied as for the data events. D0 obtains a good description of beam and 
instrumental backgrounds when overlaying events recorded with a zero bias (i.e. random) 
trigger in data. 

2'2. Single Top Production. - The Tevatron collaborations have reported evidence for 
electroweak single top quark production [2,3]. At y^s = 1.96 TeV single top quarks are 
predominantly produced via the i-channel process qb — > q't and the s-channel process 
qq' — > tb. Theoretical predictions of the production cross section are available at next- 
to-leading order (NLO), see e.g. ref. [8]. 

The LO kinematics of the s-channel process are unchanged by NLO corrections; there- 
fore it is sufficient to employ a LO MC generator to simulate the s-channel process and 
to scale the obtained cross section to the NLO expectation. On the other hand there 
are important corrections to the t-channel kinematics from the 2 — > 3 process qg — ► q'tb, 
where the gluon in the initial state produces a bb pair. The soft and collinear regime of 
this process is well modeled by the 2 — > 2 process available in LO MC with PS, using b 
quark PDFs; however, for large transverse momenta of the final state b quark, Pt(&), the 
2^3 process needs to be considered explicitly. 

Both CDF and DO use a procedure in which the 2^2 and the 2 — > 3 processes are 
generated separately and the phase space overlap between the two processes is removed 
by hand [9]. The NLO cross section ctnlo as a function of Pr(b) is then constructed as 
follows: 

(1) ^NLO = K ■ cr 2 ^2,PYTHiA| pr ( 6 ) <p o, + 0-2^3| PT (6)>pO , 

which is illustrated in fig. Q] and implies the following steps: 

• The 2 — > 2 and the 2^3 processes are generated separately, utilizing the MadEvent 
generator [12] in the case of CDF, and the SingleTop generator [9] at D0. 

• The 2 — ► 2 process is scaled up by a factor of K such that the cross section for 
both processes matches the NLO cross section. 

• The soft and the hard regime are separated by a simple cut on pr{b) at pj>, i.e. 
all events with pt{V) < p^p are taken from the 2 — > 2 process and all events with 
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Fig. 1. - Constructing the top t-channel cross section as a function of the b quark transverse 
momentum from the 2-^2 and the 2-^3 process [9-11]. 



Pr(b) > Pt are taken from the 2^3 process. The value of p 1 ^ is chosen to ensure 
a smooth transition of the differential cross section as a function of pr{b) between 
the soft and the hard regime (around 20 GeV/c) [10,11]. 

• The full event kinematics are then validated against NLO predictions from the 
ZTOP code [8] . As shown in fig. [2 both collaborations find good agreement between 
the MC samples and the NLO prediction. 

The above procedure represent a "pragmatic" solution that is sufficient for the current 
precision of the single top measurements. In the future it is desirable to treat t-channel 
single top production consistently in a full NLO MC simulation, e.g. in the framework 
of MC@NLO [13]. 



3. Monte Carlo Simulations and Systematic Uncertainties 

CDF and D0 have established their (separate) ways of assigning systematic uncer- 
tainties to top physics results over the course of Tevatron Runs I and II. A significant 
fraction of the systematic uncertainties is related to the MC simulation. For example, 
for the most recent Tevatron top mass combination, out of the total relative uncertainty 
of 0.8%, 0.3% can be attributed to MC related effects [1]. In the following we will review 
the standard treatment of MC related systematic uncertainties by CDF and DO and 
comment on recent progress on this subject. 

3T. Signal Model. - The two main all-purpose LO MC generators with PS, Pythia [4] 
and HERWIG [14], have different hadronization models and different tunings for the 
underlying event (i.e. interactions of partons inside the colliding proton and antiproton 
other than the ones involved in the hard scattering process). The treatment of initial 
state radiation (ISR) and final state radiation (FSR) in Pythia and HERWIG is similar, 
but not identical. 

CDF extracts MC model uncertainties for many observables from a comparison of 
Pythia and HERWIG. For the first years of Tevatron Run II, the DO collaboration 
did not employ the HERWIG MC generator and therefore did not compare HERWIG to 
Pythia. Instead, DO evaluates uncertainties due to the b quark fragmentation model, 
a major source of MC model uncertainty in b rich top events. 

CDF and DO also follow different approaches to evaluating uncertainties due to ISR 
and FSR. CDF estimates the uncertainty from special MC samples in which the Pythia 
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Fig. 2. - Comparison of top quark transverse momentum distributions in i-channel single top 
production. Left: CDF comparison of MadGraph and ZTOP. Right: D0 comparison of SingleTop 
and ZTOP (after ref. [9]). 



parameters controlling the amount of ISR and FSR ( Aqcd and Q 2 scales) are varied. This 
procedure is sensitive to differences in the soft radiation. D0 extracts the systematic 
uncertainties due to ISR and FSR from reweighting the jet multiplicity spectrum in 
the MC simulation (ALPGEN ti+ 0-2 partons) in a control region to the multiplicity 
spectrum observed in data. With this procedure, D0 mainly probes differences in hard 
radiation. 

Both collaboration evaluate uncertainties due to the PDF choice in similar ways. For 
any given observable, e.g. the acceptance for ti events, the uncertainty is obtained from 
two main sources, the change in the observable when varying eigenvectors in the space 
of the PDF fit parameters and the difference in the observable when calculated using 
CTEQ [7] versus MRST [15] PDF sets. 



3'2. Jet Energy Scale. - In many top analyses, uncertainties in the calibration of 
the jet energy scale (JES) are among the leading systematic uncertainties. The CDF 
procedure for JES corrections is detailed in ref. [16]. The JES correction is obtained e.g. 
from dijet and photon-jet balance, and the uncertainty of the correction is estimated, 
among others, from a comparison of Pythia and HERWIG. Note that it is important 
to avoid double counting of MC uncertainties in the signal MC model and the JES, as 
they are both derived in part from Pythia-HERWIG comparisons. 

Top quark decays are a significant source of b quarks. In general, jets containing 
b quarks have a different JES than light parton jets. The uncertainty of the b JES 
correction is larger than for light jets, because there are additional uncertainties due to 
fragmentation and color flow, and the small b quark jet sample sizes make a calibration 
in data less precise. 

In recent years the JES uncertainty has been reduced significantly in top mass analyses 
by employing an in situ JES calibration technique [17]. In the tt —> Ivbqqb decay 
channel, the mass of the top quarks and the mass of the hadronically decaying W boson 
are extracted simultaneously. The comparison of the measured W mass with the world 
average W mass is used to calibrate the JES. In this case, the JES uncertainty reduces 
to a residual uncertainty that covers the dependence on the jet kinematics. 
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Fig. 3. - Comparison of the ii transverse momenta in 1.9 fb x of data with different MC gener- 
ators. Top left: Pythia. Top right: HERWIG + JIMMY. Bottom left: ALPGEN + Pythia (tt 
+ jets). Bottom right: MCQNLO. 



3"3. Recent Developments. - The Tevatron collaborations are in the process of revis- 
iting their procedures for assigning systematic uncertainties. This becomes important 
especially for precision measurements of the top quark mass. Starting from two joint 
workshops in 2007, CDF and D0 are working towards a better understanding of each 
other's procedures and a common approach for assigning systematic uncertainties. 

One example of a long-standing disagreement between Tevatron data and MC simula- 
tions lies in the transverse momentum distribution of tt pairs. Fig. [3] shows a comparison 
of CDF data with different MC generators. None of the generators shows good agreement 
with the data. The HERWIG-based generators (HERWIG + JIMMY [18], MC@NLO) 
show better agreement with the data than the PYTHiA-based generators (Pythia, ALP- 
GEN + Pythia). Generators that include higher orders (ALPGEN + Pythia, MC@NLO) 
show the same disagreement as LO generators. These findings indicate that the source 
of the disagreement lies in the PS part of the MC generators. Note that the influence of 
the disagreement on the top quark mass uncertainty is small. 



4. Monte Carlo Simulations for Special Analyses 

In many analyses of top quark properties, the standard MC tools for tt production do 
not provide all the features required for the analysis. This is especially true for CDF's 
main generator Pythia v6.2. In this section we will discuss examples of Tevatron top 
properties analyses with special MC requirements. 
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4T. W Helicity in Top Quark Decays. - Measurements of the helicity of W bosons in 
top quark decays are tests of the V— A structure of the tWb vertex. The recent Tevatron 
analyses [19, 20] measure the W helicity via the distribution of the angle cos 9* between 
the charged lepton (or down-type quark) and the top boost direction in the W rest frame. 

To compare the acceptances for ti processes with different W helicities, the cos 9* 
distributions obtained from ALPGEN + Pythia (D0) or Pythia (CDF) are reweighted 
according to helicities other than the standard model (SM) prediction. For linearity 
checks of the W helicity measurement, CDF has also used the MadEvent generator, 
which changes the underlying left-handed and right-handed couplings rather than the 
direct observable cos9* . Systematic uncertainties of the W helicity measurement are 
obtained from comparing the SM helicities predicted by different MC generators, e.g. 
Pythia versus HERWIG in the case of CDF. 

4"2. Flavor Changing Neutral Currents in the Top Sector. - Flavor changing neutral 
current (FCNC) interactions of top quarks are strongly suppressed in the SM, so that 
any FCNC signal would be an indication for physics beyond the SM. CDF has searched 
for the FCNC t — ► Zq [21] in ti decays, and D0 has searched for single top production 
via the FCNC q -> tg [22]. 

The collaborations have chosen different approaches to determining the FCNC signal 
acceptances. CDF obtains the t — > Zq decay from Pythia and reweights the resulting 
isotropic cos 9* distribution (defined analogously to sect. I4'ip according to the expec- 
tations. This approach stays close to experimental observables, which can also be seen 
by the result, which is a limit on the branching fraction for the t — > Zq decay. The D0 
analysis utilizes the CompHEP MC generator [23] to modify the top quark couplings. 
The D0 approach is closer to theoretical calculations, and consequently, the result is 
quoted as a limit on the qtg coupling. 

4'3. Fraction of ti Events Produced in Gluon-Gluon Fusion. - NLO QCD calculations 
predict that at the Tevatron, 85% of all ti pairs are produced via qq annihilation and 
15% are produced via gg fusion. CDF has recently published first experimental tests of 
this prediction [24] . One of the CDF analyses measures the fraction of ti produced via 
gg fusion by training an artificial neural network on the kinematics of ti production and 
decay. While the main sensitivity of this method comes from the production angles and 
velocities of the top quarks, the dependence on ti spin correlations, which is encoded in 
the decay kinematics, also contributes to the sensitivity. 

The standard CDF MC tool Pythia is missing two important ingredients for this 
analysis. Pythia as a LO MC generator only generates 5% gg — ► ti, and it does not 
include ti spin correlations. Therefore the analysis makes use of the HERWIG generator, 
also using a CDF-internal extension of HERWIG named GGWIG that allows the user to 
adjust the relative fraction of gg — ► ti. The systematic uncertainty of the measurement 
due to the difference between HERWIG and Pythia is small; however, the systematic 
difference between HERWIG and the NLO generator MC@NLO results in a 20% loss in 
sensitivity, see fig. [4] This is one of the largest systematic uncertainties and indicates the 
importance of NLO effects for the analysis. 

4'4. Charge (Forward Backward) Asymmetry. - Both Tevatron collaboration have re- 
cently published measurements of the charge asymmetry in ti production (which trans- 
lates into a forward-backward asymmetry if one assumes CPT invariance) [25,26]. The 
asymmetry vanishes at LO, and NLO calculations predict an asymmetry of 4-5% [27], 
which is confirmed by the MC@NLO generator, which shows a 3.8% asymmetry. 
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Fig. 4. - Feldman-Cousins (FC) bands at 68% confidence level for the measurement of the frac- 
tion of ti production via gg fusion. The FC band for statistical uncertainties plus uncertainties 
due to the comparison of HERWIG with MC@NLO is compared with the FC band for statis- 
tical uncertainties only. At the fitted gg fraction of —7.5%, the limit on the true gg fraction is 
deteriorated by approximately 20%. 



The net asymmetry obtained from LO MC generators is zero; however, there may 
be asymmetries generated from color flow for specific kinematic configuration. This is 
illustrated in fig. [5] which shows a comparison of the asymmetry as predicted in Pythia, 
ALPGEN, and MC@NLO as a function of the rapidity difference between the t and the 
t quark. The asymmetries are also subject to NLO corrections. For example, for ti + jet 
production, a LO asymmetry of —8% is reduced to —1.5% [28]. In summary, there is 
currently no single MC generator that covers all features of the charge (FB) asymmetry. 



5. Summary and Outlook 

The CDF and D0 collaborations have entered an era of precision top quark physics. 
While more and more sophisticated analysis techniques are being employed, both collab- 
orations favor a "conservative" and "pragmatic" approach towards MC simulations: 

• Prefer well validated MC tools (mainly based on Pythia) over "cutting edge" MC 
technology. 

• Prefer data-driven methods whenever available. 

• Adapt and re-use well understood MC samples, e.g. by applying reweighting tech- 
niques. 

In addition there are often technical reasons for the above approach. Integrating new 
MC generators into the existing CDF and D0 frameworks is time intensive and error 
prone, despite close collaboration with the authors of the MC codes. With the current 
computing and person power available at the Tevatron, the cycle of software integration, 
event generation, validation, and tuning of a new MC generator can easily take many 
months. 

The Tevatron collaborations are actively working on common approaches to top MC 
generators. The final goal is a common treatment of MC-related systematic uncertainties, 
primarily for the top quark mass measurements, but also for other top analyses. 
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Fig. 5. - Charge (forward-backward) asymmetry as a function of the rapidity difference between 
the t and the i quark for the Pythia, ALPGEN, and MC@NLO generators. 
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